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Observations suggest that 7 -ray bursts (GRBs) are produced by the dissipation of the kinetic energy 
of a relativistic fireball. We show that a large fraction, > 10%, of the fireball energy is expected to 
be converted by photo-meson production to a burst of ~ 10^'*eV neutrinos. A km^ neutrino detector 
would observe at least several tens of events per year correlated with GRBs, and test for neutrino 
properties (e.g. flavor oscillations, for which upward moving r’s would be a unique signature, and 
coupling to gravity) with an accuracy many orders of magnitude better than is currently possible. 
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Recent observations of 7 -ray bursts (GRBs) suggest 
that they originate from cosmological sources 0 (see, 
however, ||]). General phenomenological considerations 
indicate that the bursts are produced by the dissipation 
of the kinetic energy of a relativistic expanding fireball 
(see for reviews). The physical conditions in the dis¬ 
sipation region imply 0 that protons may be Fermi ac¬ 
celerated in this region to energies > 10^°eV. Further¬ 
more, the spectrum and flux of ultra-high energy cosmic 
rays (above 10^®eV) are consistent with those expected 
from Fermi acceleration of protons in cosmological GRBs 
1^. We show in this Letter that a natural consequence 
of the dissipative fireball model of GRBs is the conver¬ 
sion of a significant fraction of the fireball energy to an 
accompanying burst of ~ lO^^eV neutrinos, created by 
photo-meson production of pions in interactions between 
the fireball 7 -rays and accelerated protons. The neu¬ 
trino burst is produced by interaction with protons with 
energies much lower than ~ 10 ^*^eV, the maximum accel¬ 
eration energy. As shown below, lO^^eV protons interact 
with the ~ IMeV photons carrying the bulk of 7 -ray 
energy to produce ~ 10 ^‘^eV neutrinos. 

The rapid variability time, ~ 1ms, observed in some 
GRBs implies that the sources are compact, with a lin¬ 
ear scale To ^ lO^cm. The high luminosity required for 
cosmological bursts, ~ lO^^erg s“^, then results in an 
optically thick (to pair creation) plasma, which expands 
and accelerates to relativistic velocities |^. The hard¬ 
ness of the observed photon spectra, which extends to > 
lOOMeV, implies that the 7 -ray emitting region must be 
moving with a Lorentz factor T of order 100 0, and con¬ 
stitutes independent evidence for ultra-relativistic out¬ 
flow. The high energy density in the source would result 
in complete thermalization and a black-body spectrum 
1^, in contrast with observations. To overcome this prob¬ 
lem, Rees and Meszaros suggested that 7 -ray emission 
results from the dissipation at large radius of the kinetic 
energy of the relativistic ejecta. Such dissipation is ex¬ 
pected to occur due to a collision with the inter-stellar 
medium |^, or due to internal collisions within the ejecta 
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Paczyhski and Xu suggested g, that 7 -rays are emit¬ 


ted by the decay of neutral pions, which are produced in 
pp collisions once the kinetic energy is dissipated through 
internal collisions. In this case, an accompanying burst 
of 30GeV neutrinos is expected due to the decay of 
charged pions. It is not clear, however, whether this 
model can account for the observed spectra. A more con¬ 
servative mechanism for 7 -ray production is the emission 
of synchrotron radiation (possibly followed by inverse- 
Gompton scattering) by relativistic electrons accelerated 
in the dissipation shocks |^,0. In this case, the pro¬ 
ton density in the wind is too low to allow significant 
conversion of energy to neutrinos through pp collisions. 

In the region where electrons are accelerated, protons 
are also expected to be shock accelerated. This is sim¬ 
ilar to what is thought to occur in supernovae remnant 
shocks, where synchrotron radiation of accelerated elec¬ 
trons is the likely source of non-thermal X-rays (recent 
ASGA observations give evidence for acceleration of elec¬ 
trons in the remnant of SN1006 to lO^^eV fl^), and 
where shock acceleration of protons is believed to pro¬ 
duce cosmic rays with energy extending to ~ lO^^eV 
(see, e.g. for review). The spectrum of ultra-high 
energy cosmic rays (above I0^®eV) is consistent [|| with 
that expected from Fermi acceleration of protons in cos¬ 
mological GRBs, and the flux is consistent with this sce¬ 
nario Q provided the efficiency with which kinetic en¬ 
ergy is converted to accelerated protons is comparable to 
the efficiency with which energy is converted to acceler¬ 
ated electrons (and hence to 7 -rays). We derive below 
the expected spectrum and flux of high energy neutri¬ 
nos, produced by photo-meson interactions between the 
wind 7 -rays and shock-accelerated protons, and discuss 
the implications for high energy neutrino astronomy. 

Neutrino production in dissipative wind models of 
GRBs. We consider a compact source producing a wind, 
characterized by an average luminosity L ~ lO^^erg s“^ 
and mass loss rate M = LjrjC^. At small radius, the wind 
bulk Lorentz factor, T, grows linearly with radius, until 
most of the wind energy is converted to kinetic energy 
and r saturates at T 77 ^ 100. Variability of the source 
on time scale At, resulting in fluctuations in the wind 
bulk Lorentz factor T on similar time scale, would lead to 
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internal shocks in the ejecta at a radius r ^ ~ F^cAt. 

We assume that internal shocks reconvert a substantial 
part of the kinetic energy to internal energy, which is then 
radiated as 7 -rays by synchrotron and inverse-Compton 
radiation of shock-accelerated electrons. 

The photon distribution in the wind rest frame is 
isotropic. Denoting by nj{e^)dej the number density of 
photons in the energy range to in the wind 

rest frame, the fractional energy loss rate of a proton with 
energy Cp in the wind rest frame due to pion production 
is 


~ r^/Fct^. The energy density in the BATSE range, 
U-y, is related to the luminosity by = 47 rr^r^cC/..y. 
Using this relation in (i, is given by 
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Here, = 10 ®^L.y, 5 ierg s T = SOOFaoo, At = 
10 “^AtmsS, and the proton break energy is 
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where Fp = Cp/mpC^, CT.n.(e) is the cross section for pion 
production for a photon with energy e in the proton rest 
frame, ^(e) is the average fraction of energy lost to the 
pion, and eo = O.lSGeV is the threshold energy. The 
GRB photon spectrum is well fitted in the BATSE range 
(30 KeV-3 MeV) by a combination of two power-laws, 
n{ej) oc eT^ with different values of (3 at low and high 
energy |]^. The break energy (where (3 changes) in the 
observer frame is typically e°^- ~ IMeV, with /3 ~ 1 at 
energies below the break and /3 ~ 2 above the break. 
Hereafter we denote quantities measured in the observer 
frame with the super-script “ob.” (e.g., e°^- = Te^b)- 

The second integral in i) may be approximated by 
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where U-y is the photon energy density (in the range cor¬ 
responding to the observed BATSE range) in the wind 
rest-frame, /3 = 1 for e < e^b and /3 = 2 for e > e^b- 
The main contribution to the first integral in ( 0 ) is from 
photon energies e ~ Cpeak = 0.3GeV, where the cross 
section peaks due to the A resonance. Approximating 
the integral by the contribution from the resonance we 
obtain 
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= 1.3 X eV. (5) 

Thus, for parameters typical of a GRB producing wind, a 
significant fraction of the energy of protons accelerated to 
energies larger than the break energy, ~ 10^®eV, would 
be lost to pion production. Note, that since the flow is 
ultra-relativistic, the results given above are independent 
of whether the wind is spherically symmetric or jet-like, 
provided the jet opening angle is > 1/F (for a jet-like 
wind, L is the luminosity that would have been produced 
by the wind if it were spherically symmetric). 

Since the constraint on the bulk Lorentz factor, F > 
100 , is derived from the requirement that the wind be 
optically thin to pair production for photons with ob¬ 
served energy ^ lOOMeV, it is useful to express as 

a function of the pair production optical depth A 

test photon with energy et may produce pairs in inter¬ 
actions with photons with energy exceeding a threshold 
eth, determined by ctCth = 2 (meC^)^/(l — cos 6 >), where 9 
is the angle between the photons propagation directions. 
Pair production interactions involve photons with energy 
much higher than the break energy, e-^b- A test photon 
with observed energy e°'^- = 100 MeV, for example, has 
an energy ^ IMeV in the wind rest frame, and there¬ 
fore interacts mainly with photons of similar energy to 
produce pairs. Assuming that the spectrum of photons 
produced in the dissipation region extends as n(e) oc 
from the BATSE range to (observed) energies well above 
lOOMeV, the mean free path for pair production (in the 
wind rest frame) for a photon of energy et is 
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Here, dpeak — 5 x 10“^®cm^ and ^peak — 0.2 are the values 
of a and ^ at e = Cpeak, and Ae ~ 0.2GeV is the peak 
width. 

The energy loss of protons due to pion production 
is small during the acceleration process Q. Once ac¬ 
celerated, the time available for proton energy loss by 
pion production is comparable to the wind expansion 
time as measured in the wind rest frame, td ~ rd/Fc. 
Thus, the fraction of energy lost by protons to pions is 


Here we have used a constant cross section, 3ctt /16 where 
(JT is the Thomson cross section, above the threshold eth 
[The cross section drops as log(e)/e for e eth; however, 
since the number density of photons drops rapidly with 
energy, (i is a good approximation]. 

The optical depth is given by = Vd/Tly^. Using 
(^, we obtain 

, , min(e°*’-, e°U) 

Uef) = 0.15r.,^(ef ■ = lOOMeV) . (7) 
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Only a small fraction of bursts show a power law spec¬ 
trum extending to > lOOMeV. Most bursts may therefore 
have = lOOMeV) larger than unity, leading to 

higher efficiency of pion production. Furthermore, vari¬ 
ability of the source over different time scales would lead 
to dissipation shocks over a range of radii 0 , where the 
optical depth may become small only at the largest dis¬ 
sipation radii. While photons can escape only from radii 
where the optical depth is small, neutrinos can escape 
from almost any depth. Thus, most of the burst energy 
may actually come out in neutrinos. For F = 100 and 
L-^ = lO^^erg s“^, for example, = lOOMeV) ~ 1 

for Td 10^‘^cm, where internal shocks result from vari¬ 
ability over 0.1s time scale. Variability on shorter time 
scales results in collisions at smaller radii with larger pion 
production efficiency. The conversion to high energy pi- 
ons of a significant fraction, > 20%, of the energy of 
protons accelerated to energy similar to or larger than 
the break energy Cpb, would be avoided only if F » 100. 

Neutrino spectrum and flux. Roughly half of the en¬ 
ergy lost by protons goes into 7r° ’s and the other half 
to 7r+ ’s. Neutrinos are produced by the decay of tt+’s, 
7 r+ ^ + e'^ +Ve + Vfj_ + v^ [the large optical depth 

for high energy y’s from 7r° decay, cf. eq. (^, would 
not allow these photons to escape the wind]. The mean 
pion energy is 20% of the energy of the proton producing 
the pion. This energy is roughly evenly distributed be¬ 
tween the 7r+ decay products. Thus, approximately half 
the energy lost by protons of energy Cp is converted to 
neutrinos with energy ~ O.OfiCp. Eq. (|^ then implies 
that the spectrum of neutrinos above Cj,}, = 0.05ep{, fol¬ 
lows the proton spectrum, and is harder (by one power 
of the energy) at lower energy. The break energy is 
= 5 X 10i^F|oo(e°)]-^jgy)-ieV. For a power law dif¬ 
ferential spectrum of accelerated protons n(ep) oc e“^, 
as typically expected for Fermi acceleration and which 
would produce the observed spectrum of ultra-high en¬ 
ergy cosmic rays [|j, the differential neutrino spectrum 
is n(ei/) oc e““ with a = 1 below the break and a = 2 
above the break. The spectrum may be modified above 
the energy where are produced by decay of muons with 
life time F^^t^ (where F^ is the muon Lorentz factor and 
tp = 2 X 10“®s its rest life time) comparable to the wind 
dynamical time rd/Tc, i.e. above e°^' = mpC^rd/^tpC « 
T'^mpC^At/3tp = 5 X 10^®r|gQAtniseV. Muons producing 
Vp with higher energy may adiabatically lose a significant 
part of their energy before decaying (Other energy loss 
processes, e.g. inverse-Compton scattering, can be shown 
to be less important). 

Most of the neutrino energy is carried by neutrinos 
with energy close to the break energy, e°)]- ~ lO^^eV. 
The neutrino flux depends on the relative efficiency with 
which the wind kinetic energy is converted to acceler¬ 
ated protons, compared to the efficiency with which en¬ 
ergy is converted to accelerated electrons (and therefore 


to 7 -rays). If cosmological GRBs are the sources of ultra- 
high energy cosmic rays, then the efficiency of converting 
energy to y-rays and to accelerated protons should be 
similar. The energy production rate required to produce 
the observed flux of ultra-high energy cosmic-rays, as¬ 
suming that the sources are cosmologically distributed, 
is ^ 4 X lO^'^'^erg Mpc“^yr“^ over the energy range 10^®- 
lO^^eV, comparable to the rate of energy production by 
GRBs as y-rays in the BATSE range ||]. For a pro¬ 
ton generation spectrum n(ep) oc e“^, similar energy is 
contained in equal logarithmic energy intervals, implying 
that the total energy converted to accelerated protons 
(i.e. over the energy range FpTOpC^ ~ lO^^-lO^^eV) is a 
few times that converted to y-ray energy in the BATSE 
range. We therefore assume below that the conversion 
efficiency is similar for electrons and protons, leading to 
energy production rate above the proton energy break, 
^ 10^®eV, E 4 X 10^‘^erg Mpc~^yr“^ {E depends only 
logarithmically on the value of the break energy). 

The present day neutrino energy density due to GRBs 
is approximately given by U^, « 0.5fTr{epb)tHE, where 
tn ~ lO^^yr is the Hubble time. The neutrino flux is 
therefore approximately given by 
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where E = 10^'^if44erg Mpc“^yr“^. The high energy 
neutrinos predicted in the dissipative wind model of 
GRBs may be observed by detecting the Cherenkov light 
emitted by high energy muons produced by neutrino in¬ 
teractions below a detector on the surface of the Earth 
(see 0 for a recent review). This technique works only 
for muons entering the detector from below, “upward 
moving muons”, due to the large background of down¬ 
ward atmospheric muons. The probability P^p that a 
neutrino would produce a high energy muon in the detec¬ 
tor is approximately given by the ratio of the high energy 
muon range to the neutrino mean free path. At the high 
energy we are considering, Pyp cx 10“®(ei//lTeV) [ p^ . 
Using (^, the expected flux of upward moving muons is 
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The rate is almost independent of e^b, due to the increase 
of Pvp with energy. 

The rate (^) is comparable to the background expected 
due to atmospheric neutrinos |14|. However, neutrino 
bursts should be easily detected above the background, 
since the neutrinos would be correlated, both in time and 
angle, with the GRB y-rays. A km^ neutrino detector 
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should detect each year ~ 10 to 100 neutrinos correlated 
with GRBs. Furthermore, nearby bright GRBs, although 
rare (~ 0.1 per year), would produce a burst of several 
neutrinos in a km^ detector. A burst at a distance of 
100 Mpc producing 0.4 x lO^^erg in ^ lO^'^eV neutrinos 
would produce ^ 3 km~^ upward muons. 

Implications. Detection of neutrinos from GRBs 
would corroborate the cosmological hreball scenario for 
GRB production (acceleration of protons to ultra-high 
energy by the processes discussed above is not possible 
if GRBs are Galactic ||^). Neutrinos from GRBs could 
be used to test the simultaneity of neutrino and photon 
arrival to an accuracy of ~ 1 s (~ 1 ms for short bursts), 
checking the assumption of special relativity that photons 
and neutrinos have the same limiting speed [The time de¬ 
lay for neutrino of energy lO^^eV with mass rrii, traveling 
100 Mpc is only ^ 10“^^(m,y/10 eV)^s]. These observa¬ 
tions would also test the weak equivalence principle, ac¬ 
cording to which photons and neutrinos should suffer the 
same time delay as they pass through a gravitational po¬ 
tential. With 1 s accuracy, a burst at 100 Mpc would re¬ 
veal a fractional difference in limiting speed of 10“^®, and 
a fractional difference in gravitational time delay of order 
10“® (considering the Galactic potential alone). Previous 
applications of these ideas to supernova 1987A (see Hi 
for review), where simultaneity could be checked only to 
an accuracy of order several hours, yielded much weaker 
upper limits: of order 10“® and 10“^ for fractional dif¬ 
ferences in the limiting speed 0 and time delay jT^ ] 
respectively. 

The model discussed above predicts the production of 
high energy muon and electron neutrinos with a 2:1 ratio. 
If vacuum neutrino oscillations occur in nature, then neu¬ 
trinos that get here should be almost equally distributed 
between flavors for which the mixing is strong. In fact, 
if the atmospheric neutrino anomaly has the explanation 
it is usually given, oscillation to i/r’s with mass ~ 0.1 eV 
0, then one should detect equal numbers of z/^’s and 
i^r’s. Upgoing r’s, rather than ^’s, would be a distinc¬ 
tive signature of such oscillations. Since z^r’s are not ex¬ 
pected to be produced in the fireball, looking for upgoing 
r’s would be an “appearance experiment” (z/^’s may be 
produced by photo-production of charmed mesons; How¬ 
ever, the high photon threshold, ~ SOGeV, and low cross- 
section, ~ 1/rb [0, for such reactions imply that the ra¬ 
tio of charmed meson to pion production is ~ 10“^). To 
allow flavor change, the difference in squared neutrino 
masses. Am?, should exceed a minimum value propor¬ 
tional to the ratio of source distance and neutrino energy 
[ 0 . A burst at 100 Mpc producing lO^'^eV neutrinos 
can test for Am^ > 10“^®eV^, 5 orders of magnitude 
more sensitive than solar neutrinos. Note, that due to 
the finite pion life time, flavor mixing would be caused 
by de-coherence, rather than by real oscillations, for neu¬ 
trinos with masses > O.leV. 

The may resonantly oscillate to z/g’s (the MSW 


effect) as they escape the fireball, provided the elec¬ 
tron number density exceeds the MSW resonance density 
(e.g. 0). The hreball electron density is approximately 
rie ~ L/Airr^T'^mpC^ ~ 10^°-10^^cm“^, implying that 
the MSW effect could be important only for neutrino 
masses smaller than ^ 10“^^eV^. For such low masses, 
the spatial scale of the hreball implies that rze changes too 
rapidly to allow the use of the adiabatic approximation, 
and that the inhuence of the MSW effect would depend 
on the detailed structure of the hreball. 
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